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It has been a long-standing goal to create magnetism in a nonmagnetic material by manipulating 
its structure at the nanometer scale. This idea may be realized in graphitic carbon: evidence 
suggests magnetic states at the edges of graphene ribbons and at grain boundaries in graphite. Such 
phenomena have long been regarded as unlikely in silicon because there is no graphitic bulk phase. 
Here we show theoretically that intrinsic magnetism indeed exists in a class of silicon surfaces whose 
step edges have a nanoscale graphitic structure. This magnetism is intimately connected to recent 
observations, including the coexistence of double- and triple-period distortions and the absence of 
edge states in photoemission. Magnetism in silicon may ultimately provide a path toward spin-based 
logic and storage at the atomic level. 



Individual atoms with an odd number of electrons ex- 
hibit a magnetic moment from the spin of the unpaired 
electron. Elements with an even number of electrons, 
such as carbon and silicon, can also reveal unpaired elec- 
trons when covalent bonds are broken. In defective or 
disordered group-IV solids the single electrons that oc- 
cupy dangling bonds can be detected in magnetic reso- 
nance experiments, but the lack of long-range structural 
order precludes a magnetically ordered state. 

Bonds are also broken at the surfaces and interfaces 
of solids, but covalent materials usually reconstruct to 
eliminate dangling bonds or doubly occupy them with 
electrons of opposite spin. Exceptions to this behav- 
ior provide fertile ground for magnetism in nonmagnetic 
materials. Two widely discussed examples are based on 
graphitic carbon: highly oriented pyrolytic graphite ex- 
hibits ferromagnetic order arising from two-dimensional 
arrays of defect spins at grain boundaries, 1 and nanoscale 
graphene ribbons are predicted to have ferromagnetically 
ordered edge states^ 

The edges of graphene ribbons are theoretically ap- 
pealing but far from being realized experimentally with 
the required atomic smoothness. Other one-dimensional 
systems such as nanowires^ and step edges provide more 
favorable conditions for long-range structural order and 
thus are promising platforms for devices requiring or- 
dered arrays of spins. Indeed, atomically smooth steps 
with structural order extending to micrometers can be 
produced by self-assembly at vicinal silicon surfaces. 
Some of these surfaces naturally form a graphitic sil- 
icon ribbon at the step edge. Here we predict, using 
density-functional theory (DFT, see Methods), that the 
electronic ground state of these silicon surfaces is mag- 
netic. Magnetism in silicon systems having no magnetic 
elements is quite unexpected. Indeed, none of these sys- 
tems has yet been investigated for magnetism experimen- 
tally, but we show below that existing data already sup- 
port it indirectly. 

Our discussion focuses on the Si(553)-Au and Si(557)- 
Au surfaces, shown in Fig. 1. These are stepped surfaces 
stabilized by a fraction of a monolayer of Au. Crystallo- 
graphically they can be considered miscut from the flat 



Si(lll) surface by angles of 12.3° and 9.4° away from 
and toward, respectively, the (001) orientation. They are 
members of a family of Au-induced vicinal silicon surfaces 
intermediate between (111) and (001) containing several 
structural motifs in common^ These commonalities sug- 
gest that magnetism in silicon may be widespread. 

The critical structural motif for magnetism in these 
systems is a nanoscale honeycomb strip of graphitic sil- 
icon, just a single hexagon wide, that forms the edge of 
the step£ This honeycomb chain was first observed in 
reconstructions of flat Si(lll) induced by submonolayer 
coverages of alkali, alkaline-earth, and some rare-earth 
metals, as well as Ag and Au£ In these flat surfaces the 
orbitals of the outer honeycomb atoms form filled bands, 
either by accepting electrons from the adsorbates or by 
forming partially covalent bonds. At a step edge a new 
possibility arises in which these orbitals are singly occu- 
pied, fully spin polarized, and magnetically ordered. 

The two surfaces studied here have step edges with 
opposite orientation. The DFT ground state (see Meth- 
ods) depicted in Fig. la for Si(553)-Au shows every third 
Si atom at the step edge to be completely spin polar- 
ized while all other atoms are negligibly polarized. On 
Si(557)-Au (Fig. lc) every second Si atom at the step is 
fully polarized as are all the terrace restatoms. In both 
systems the interaction between neighboring spins along 
a step is antiferromagnetic. 

These predictions depend on an accurate knowledge 
of the atomic structure of each surface. The model 
of Si(557)-Au shown in Fig. lc is well-established from 
DFT calculations^ and X-ray diffraction^ The energy 
bands measured in Si(557)-Au agree well with the calcu- 
lated bands arising from the Au chains. A small band 
splitting^— originates from the spin-orbit interaction of 
the Au atomsi^i^ In this "Rashba effect" — which is un- 
related to the magnetism we predict here — the two split 
bands have a spin polarization of 100% at any given mo- 
mentum in reciprocal space, but this vanishes when in- 
tegrated over all momenta and leaves no net spin polar- 
ization in real space. 

Existing models for the structure of the Si(553)- 
Au surface do not reproduce the experimental band 
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FIG. 1: Ground state structure and lowest energy spin configuration of two magnetic silicon surfaces, a— b, 

Si(553)-Au in its antiferromagnetic ground state. Yellow atoms are Au, all others are Si. Each terrace contains a Au double 
row and a graphitic Si honeycomb chain (green) at the step edge. Every third Si (red, blue) at the step has a spin magnetic 
moment of one Bohr magneton (S = 1/2, arrows) from the complete polarization of the electron occupying the dangling-bond 
orbital. The sign of the polarization alternates along the step. The periodicity along the step is tripled by a small downward 
displacement of the magnetic atoms. The periodicity along the two Au rows is doubled, with alternating short and long bonds 
Oi. The full 1x6 unit cell is outlined, c, Si(557)-Au in its antiferromagnetic ground state. The step direction of Si(557) is 
opposite to that of Si(553) as shown by the schematic outlines. Nevertheless this surface is also magnetic: every second Si (red, 
blue) at the step has a magnetic moment of one Bohr magneton, as do all of the Si restatoms (blue, red) on the terrace. 



structure . 14 ' 15 We demonstrate below that the model pro- 
posed in Fig. la correctly predicts not only the band 
structure, but also the previously unexplained coexis- 
tence of two periodicities observed in scanning tunneling 
microscopy (STM). We also show that these predictions 
are inextricably linked with magnetic order. Their exper- 
imental confirmation constitutes indirect but compelling 
evidence for the existence of magnetism as well. 

The Si(553)-Au surface was first observed in Ref. [l6| 
and its unusual one-dimensional electronic properties 
have since been extensively studied. Two experimen- 
tal findings are especially striking. (1) STM images 
acquired at room temperature show alternating bright 
and dim rows, each with the periodicity of the Si sur- 
face lattice constant a ^r— At temperatures below 50 
K these rows separately develop higher-order periodic- 
ity: a tripled period (3ao) for the bright rows and a dou- 
bled period (2ao) for the dim rows. These effects have 
variously been attributed to periodic lattice distortions^ 
and to charge-density waves^ but their true microscopic 
origin has remained elusive. (2) Angle-resolved pho- 
tocmission spectra (ARPES) obtained at 160-220 K re- 
veal three metallic bands, all with their minima at the 
zone boundary ZB|| of the lxl Brillouin zoneJ^ The two 
lower bands are slightly split by a small momentum split- 
ting 5k F = 0.05 A" 1 ^ and cross the Fermi level about 
halfway to the zone center. Detailed analysis of this 
crossing shows that the splitting is due to the spin-orbit 



interaction.— The upper band crosses the Fermi level 
about one-fourth of the way to the zone center without 
a detectable splitting. 

All of these observations are correctly predicted by the 
proposed model of Si(553)-Au. Figure 2 shows the theo- 
retical simulated STM image for the model in Fig. la in 
its antiferromagnetic ground state. The bright row arises 
from Si atoms at the step edge and indeed has tripled pe- 
riodicity 3ao. The highest peaks within this row are from 
the spin-polarized atoms. This appears counterintuitive 
since these atoms are actually 0.3 A below the height of 
their non-polarized neighbors. Moreover, the dangling- 
bond orbital on this atom is occupied by a single electron 
rather than two. But the average energy of the occupied 
spin-polarized orbital is energetically closer to the Fermi 
level and as a result the topography appears higher in 
filled-state tunneling. 

The dim second row visible in Fig. 2a arises from 
the right-hand leg of the ladder formed by two rows 
of Au atoms. Its doubled periodicity 2ao is due to the 
dimerization created when the rungs of this ladder rotate 
with alternating signs along a row. A similar dimeriza- 
tion occurs in a closely related system, Si(lll)-(5x2)- 
Au, where it it is driven by electron doping from Si 
adatoms and the consequent formation of a surface band 
gapi^i For Si(lll)-(5x2)-Au the dimerization parameter 
d = (ai — ao)/ao, where ao is the surface lattice constant 
and a\ is defined in Fig. lb, has the value d = 0.14, much 
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FIG. 2: Scanning tunneling microscopy of Si(553)- 

Au. a, Theoretical STM topography showing the tip height 
at constant current for tunneling into empty surface states at 
bias voltage +0.5 V. Inset: Electron spin density 1 A above 
the step edge. Red and blue contours represent positive and 
negative polarization, indicating antiferromagnetic alignment 
of the red and blue Si atoms the step edge. Contours are 
separated by a factor of two. b, Line scans at the locations 
marked in a, showing the coexisting 3ao and 2ao periodicity 
along the Si step edge and Au chain, respectively. 



larger than the value d = 0.04 for Si(553)-Au. But even 
this small dimerization is sufficient to produce a doubled 
periodicity that is well-resolved in the theoretical STM 
image. 

The spin polarization of the electron density, shown in 
the inset to Fig. 2a, is strongly localized at every third 
Si step-edge atom. Integrating this spin density shows 
that the red and blue step-edge atoms each have a spin 
moment of 1 Bohr magneton, equivalent to a fully po- 
larized electron with spin S — 1/2. The direction of the 
spin alternates along the step in this antiferromagnetic 
ground state. Moments on other atoms along the step 
are smaller by an order of magnitude. 

The connection between magnetism and the periodic- 
ity of the rows is simple and direct. If the spin polariza- 



tion is constrained to be zero then the DFT equilibrium 
geometry of the model changes: both the tripled period- 
icity along the Si rows and the doubled periodicity along 
the Au rows are completely eliminated. In this sense 
the experimental observation of triple and double peri- 
odicities constitutes strong evidence for a spin-polarized 
ground state. 

Magnetism also provides the key to understanding 
a long-standing puzzle presented by ARPES data for 
Si(553)-Au. The orbitals at the step edges of a vicinal 
surface would normally give rise to sharp electronic states 
near the Fermi level, as is indeed the case for previously 
proposed models of Si(553)-Aui 14 ' 15 But these are not ob- 
served in the ARPES data. We demonstrate below that 
spin polarization at the step edge provides the resolution 
to this puzzle. 

Figure 3a shows the theoretical bands for Si(553)-Au 
unfolded into the extended Brillouin zone, where they 
can be easily compared to those obtained from ARPES 
spectra. Since the theoretical ground state is antiferro- 
magnetic, one need only analyze the bands from a single 
spin channel. The symbols and curves identify, using 
the colors from Fig. 1, the atoms that dominate the cor- 
responding wavefunctions. States arising from the mag- 
netic atoms experience an exchange splitting of about 0.5 
eV. As a result, for spin-up states the red step-edge band 
is singly occupied and the blue step-edge band is empty. 
For spin-down states this situation is reversed, as seen in 
the density of states, Fig. 3b. 

While these magnetic step-edge states are responsi- 
ble for the 1x6 structural and magnetic periodicity of 
Si(553)-Au, the remaining features of the band structure 
have simple lxl nonmagnetic character. This is evident 
from comparing the 1x6 bands to those in Fig. 3d for 
a simpler lxl structure identical to our model but con- 
strained to unit periodicity ao and thus having no spin 
polarization. In this lxl structure there is just a single 
type of step atom and therefore only one corresponding 
band (green). The remnant of this nonmagnetic band 
is also visible in the full 1x6 bands, although there its 
spectral weight at a given momentum is spread out by 
partial hybridization with its magnetic counterparts. 

Both the 1x6 and lxl structures give rise to several 
metallic bands (yellow) with their minima at the zone 
boundary ZBm. These are the bands detected in ARPES. 
In the 1x6 case there are a total of four bands in the two 
spin channels, but the spin-up and spin-down bands are 
nearly degenerate. In the nonmagnetic lxl case there 
are also four bands but here the degeneracy is broken 
by the spin-orbit interaction, — which for simplicity was 
not included in the 1x6 calculation.— The resulting mo- 
mentum splitting 5k p = 0.04 A -1 of the lower band is 
in excellent agreement with the experimental splitting. 
The splitting of the upper band is less than half this 
value and consequently is not resolved in the ARPES 
data. The predicted Fermi-level crossings of these bands 
are in qualitative agreement with the measured values: 
the lower doublet crosses the Fermi level halfway to the 
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FIG. 3: Electronic band structure and density of states for Si(553)-Au. a— b, Bands and DOS for the ground-state 
antiferromagnetic structure of Fig. la. The full 1x6 bands are shown unfolded and parallel to the chain direction (ZBii is the 
boundary of the lxl Brillouin zone) along the line halfway to the orthogonal zone boundary ZBx . Red and blue diamonds are 
spin-polarized states from red and blue Si step-edge atoms in Fig. la. Only spin-up states are shown here — antiferromagnetic 
order implies the spin-down states are equivalent but with red and blue reversed, as seen in the DOS. Green diamonds are 
non-polarized states from green Si step-edge atoms. Yellow circles are bonding (lower) and antibonding (upper) states from 
the Au chains. Curves show schematic energy bands without hybridization, c— d, Bands and DOS for the same model when it 
is constrained to lxl periodicity, for which the ground state is not spin-polarized. Spin-orbit coupling, included in the bands 
but not in the DOS, splits each of the two Au chain states but not the Si step-edge state. 



zone center, and the upper doublet about one-fifth of the 
way. 

One feature of the nonmagnetic lxl band structure 
is not observed in the ARPES data: the green step- 
edge band that just crosses the Fermi level in Fig. 3d. 
Magnetism provides a simple explanation for the appar- 
ent absence of this band. Above a temperature scale 
set by the energy gained from spin polarization, thermal 
fluctuations between polarized and unpolarized electron 
states will broaden the step-edge bands. The magnitude 
of the broadening is roughly the exchange splitting, 0.5 
cV. The resulting reduction in intensity makes the step- 
edge bands essentially invisible compared to the sharp Au 
bands, which are not affected by the spin polarization at 
the step edge. In this sense the experimental absence of a 
sharp step-edge band constitutes additional evidence for 
a spin-polarized ground state. 

We now briefly summarize our predictions for mag- 
netism in Si(557)-Au, which are qualitatively similar to 
those for Si(553)-Au, and defer a detailed discussion. De- 
spite having a more complicated arrangement of spins, 
the surface has a simpler 1x2 periodicity which persists 
even if spin polarization is suppressed. Magnetism leads 
to quantitative but not qualitative structural changes: 
a smaller downward displacement of the magnetic step- 
edge atoms and a smaller upward displacement of the 
magnetic restatoms. Therefore we do not predict any 
higher-order periodicity to mark the onset of magnetism 
at low temperature. 

Earlier theoretical studies of Si(557)-Au included spin- 
orbit coupling but not spin polarization. 7 The result- 
ing prediction of momentum-split Au bands were in 
good agreement with experiment. But these studies also 



predicted Si step-edge and restatom bands not seen in 
ARPES i 4 ' 10 ' 12 This discrepancy is resolved in the mag- 
netic ground state: the exchange splitting pushes the two 
occupied spin bands down to 0.4-0.5 eV below the Fermi 
level, and the two unoccupied spin bands up to 0.1-0.2 
eV above the Fermi level. This exchange splitting is sim- 
ilar to the splitting in Si(553)-Au. Thermal fluctuations 
between polarized and unpolarized states will likewise 
broaden these step-edge and restatom bands, rendering 
them invisible to ARPES while leaving the momentum- 
split Au bands unaffected. 

We anticipate that spin-polarized states in silicon are 
not limited to Si(553)-Au and Si(557)-Au: the key struc- 
tural element — a graphitic step edge — is known to exist 
on other stepped silicon surfaces as welli^ Direct experi- 
mental tests of our predictions using, for example, spin- 
polarized ST M 24 i 25 or spin-polarized photoemission 2 ^ 
will be an important next step, and will provide impor- 
tant information about the nature of the magnetic order- 
ing. 

To set the stage for such tests it is helpful to esti- 
mate the magnitude of the spin interactions within a 
simple nearest-neighbor Heisenberg Hamiltonian. From 
DFT total-energy calculations of different spin configu- 
rations on Si(553)-Au we find antiferromagnetic coupling 
( J|| = 15 meV) along the steps and weaker, ferromagnetic 
coupling (J± = —0.3 meV) across the steps. We also 
find that the magnitude, and even the sign, of these cou- 
plings can be changed by doping electron or holes into the 
Si(553)-Au surface states. It is well-established that sur- 
face states on the closely related Si(lll)-Au surface can 
be electron-doped by adsorbates (e.g. silicon adatoms) 
on the surface, and that the concentration of this ad- 
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sorbate population can be controlled to some extent^ 
No such studies have yet been reported for Si(553)-Au 
or Si(557)-Au, but the possibility of tuning surface mag- 
netism using surface chemistry suggests a wealth of new 
research possibilities. 

Linear chains of spin-polarized atoms provide atom- 
ically perfect templates for the ultimate memory and 
logic, in which a single spin represents a bit. Here the 
spins are locked into a self-assembled rigid lattice on a 
substrate that is compatible with silicon technology. One 
potential application is the spin shift register recently 
proposed theoretically by Mahan. 27 This device requires 
a one-dimensional chain of identical atoms each with spin 
S = 1/2. When an additional electron is conducted along 
the chain, each spin state is shifted by one atom — the 
spin analogue of a standard shift register memory de- 
vice. The arrangment of spins is arbitrary and there is 
no requirement of long-range order. However, two other 
important criteria must be satisfied: in the ground state 
each atom must have one electron, and the atoms must 
have correlated electronic states. The first criterion is 
indeed satisfied by the spin-polarized atoms on Si(553)- 
Au and Si(557)-Au. In previous work we showed that 
the electronic states of singly-occupied silicon dangling 
bonds are indeed strongly correlated^ 

Another potential application is the storage of infor- 
mation in single magnetic atoms. 2 - Toward this end, Hir- 
jibehedin et al. recently used spin excitation spectroscopy 
to demonstrate that the orientation of an individual spin 
can be measured with an STM^S For information stor- 
age the stability of this orientation, for example from 
coupling of the spin to the substrate, is critical. This 
coupling is most likely to be substantial when the sub- 
strate is anisotropic at the atomic scale. In the experi- 
ments of Hirjibehedin this anisotropy was engineered by 
creating CU2N islands that break the fourfold symmetry 
of the Cu(OOl) substrate. In the Si(553)-Au and related 
vicinal silicon surfaces this anisotropy arises naturally at 
the step edge. The investigation of spin-lattice coupling 
in these magnetic silicon surfaces is an important next 
step toward the ultimate goal of spintronics with single 
spins. 



Methods 



First-principles total-energy calculations were used to de- 
termine equilibrium geometries and relative energies of 
the models in Fig. 1 and their variants. The calculations 
were performed in a slab geometry with six or more lay- 
ers of Si plus the reconstructed top surface layer and 
a vacuum region of at least 10 A. All atomic positions 
were relaxed, except the bottom Si layer and its pas- 
sivating hydrogen layer, until the largest force compo- 
nent on every atom was below 0.02 eV/A. Total energies 
and forces were calculated within the PBE generalized- 
gradient approximation 3 ^ to density-functional theory 
(DFT) using projector-augmented-wave potentials, as 
implemented in VAS P 32 i 33 . The plane-wave cutoff for all 
calculations was 250 eV. 

The sampling of the surface Brillouin zone was cho- 
sen according to the size of the surface unit cell and the 
relevant precision requirements. For the 1x6 reconstruc- 
tion of Si(553)-Au shown in Fig. l(a,b), the equilibrium 
geometries and total energies of nonmagnetic, ferromag- 
netic, and antiferromagnetic configurations were calcu- 
lated using 6x6 zone sampling, with convergence checks 
using 10x10 sampling. For the 1x2 reconstruction of 
Si(557)-Au shown in Fig. lc, the equilibrium geometries 
and total energies of nonmagnetic, ferromagnetic, and 
antiferromagnetic configurations were calculated using 
2x8 zone sampling, with convergence checks using 2x16 
sampling. 

Exchange coupling constants J|| and J± for Si(553)- 
Au were extracted from DFT total energies. Supercells 
containing two spins on each of two adjacent steps were 
constructed with the four possible unique spin configu- 
rations. The DFT total energies were fit to a nearest- 
neighbor classical Heisenberg Hamiltonian. The numer- 
ical reliability of the coupling constants was analyzed 
with respect to Brillouin-zone sampling, plane- wave cut- 
off, and relativistic treatment, and the dependence of 
the coupling constants upon charge doping of the surface 
states was investigated (see Supplementary Information) . 

The simulated STM images in Fig. 2 were calculated 
using the method of Tersoff and Hamanrt 3 -- For this 
filled-state image we integrated the local density of states 
(LDOS) over a 0.5-eV energy window of occupied states 
up to the Fermi level. The simulated STM topography 
under constant-current conditions was obtained by plot- 
ting the height at which the integrated LDOS is constant. 
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